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ABSTRACT 

We have searched a sample of 151 young, energetic pulsars for periodic variation in pulse time-of- 
arrival arising from the influence of planetary companions. We are sensitive to objects with masses two 
orders of magnitude lower than those detectable with optical transit timing, but we find no compelling 
evidence for pulsar planets. For the older pulsars most likely to host planets, we can rule out Mercury 
analogues in one third of our sample and planets with masses >0.4 M 0 and periods Pb <lyr in all 
but 5% of such systems. If pulsar planets form primarily from supernova fallback disks, these limits 
imply that such disks do not form, are confined to <0.1 AU radii, are disrupted, or form planets more 
slowly (>2Myr) than their protoplanetary counterparts. 

Subject headings: pulsars:general,planets and satellites:formation 


1. INTRODUCTION 


Somewhat surprisingly, the first extra-solar planets 
were detected i n orbit around the milli second pulsar 
PSR B1257-I-12 (Wolszczan & Frail 1992). Since then, 
optical observers have honed their techniques to the 
point where thousands of planets have now been detected 


around virtually all stellar classes (e.g. Wright et al. 2011 


Rowe et al. 2014). Meanwhile, i n spite ot potentially high 


sensitivity to low-mass planets (Thorsett & Phillips 1992_ 
Cordes|1993 ), the radio pulsar community has turned up 


only one further case, a super-Jupiter in a tr iple system 


within a g lobular cluster (PSR B1620—26, Sigurdsson 
eral][2003| . 


Why IS this so? First, a Galactic field pulsar planet 
must either have survived a supernova or have formed 


from the debris of the explosion. Phinney & Hansen 
(1993) give a delightful review of these “Salamander” and 
“Meninonides” scenarios, of which w e consider planet for¬ 


mation in a supe rnova fallback disk (Lin, Woosley & Bo- 
denheimer 1991| the most relevant mechanism. Here, 


reverse shock waves generated at density transitions in 
the stellar core and envelope decelera te material tha t 
falls back toward the compact remnant (Chevalier[|1989). 
The total mass, accretion rate, and intrinsic angular mo- 
mentum depend sensitively on both the initial stellar 
conditions and mixing during the explosion, but ulti¬ 
mately ^O.IMq of material m ay circularize into a rel¬ 


atively long- li ved “debris” disk (Menou, Perna fc H ern- 
quist 2001b). Wang, Chakrabarty fc Kaplan ( |2006[ ) ob¬ 


served a distinct infrared component from the anomalous 
X-ray pulsar 4U 0141-1-61, consistent with the reprocess¬ 
ing of X-ray emission by a warm disk, but other direct 
disk searches have failed to locate further examples (e.g. 
Wang et al.||2014 ). 


Secondly, pulsars are far from perfect clocks. A com¬ 
panion induces a neutron star reflex motion detectable as 
a periodic delay in the pulse time-of-arrival (TOA ). But 


pulsars suffer large amplitude, red timing noise (Shan¬ 


non & Cordes 2010) and glitches which can mask such 
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Figure 1. The spin-down power {E = 10^®P/P^ergs and 
characteristic ages (tc = P/2Ps) of the 151 pulsars in the sam¬ 
ple. The sample divisions (see main text) in Tc are indicated by 
the vertical lines. Pulsars suffering large glitches (Sh/ju > 10“^) 
during our observation span are indicated by coloured, distinct 
symbols (legend). Lines of constant magnetic field {B = 3.2 X 
10l9(pp)l/2 Q) 

are shown as dashed salmon lines. 


periodic signals. In addition, correlations between pulse 
shape and timing residuals indicate quasi-periodic state 
switching (Lyne et al. 2010 Hobbs, Lyne & Kramer 2010| 
with time-scales ot months to years, mimicking the signa- 
ture of planets. Finally, planets with very short or very 
long periods can evade detection due to finite observing 
cadences and data spans, respectively. 

To further investigate the paucity of pulsar planets, 
we have searched for periodic signals in a sample of 151 


energetic (spin-down luminosity E > 10^^ ergs pul¬ 
sars timed with the Parkes telescope in support of the 
Fermi mission. The E and characteristic age of the 
sample are shown in Figure These young pulsars tend 
to glitch regularly and are also subject to strong tim¬ 
ing noise, complicating the search for planets with pe¬ 
riods much longer than a year. However, given the an- 
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density as 


Figure 2. The deviation of TOAs of five pulsars to a simple spin- 
down model including i) to reduce low-frequency timing noise. The 
pulsars are chosen to illustrate the range of red and white noise in 
the sample, and residuals have been scaled by the indicated amount 
to facilitate comparison. 

gular momentum budget (^ 10 "^® ergs) of a fallback disk, 
any Earth-mass planets should reside within ^2AU of 
the neutron star, making our insensitivity to long peri¬ 
ods unimportant for searching for planets formed in such 
disks. Moreover, the age range of our sample precisely 
spans the lifespan of an active disk and subsequent planet 
formation. 

Below, we present a new technique for jointly mod¬ 
elling stochastic timing noise with other timing parame¬ 
ters which improves the robustness of detections of and 
UOTer limits on periodic signals. We derive upper limits 
( 1 ^ for a range of periods for each member of our sam¬ 
ple and present the combined constraints on the young 
pulsar population. We interpret the limits and their im¬ 
plications for supernova fallback debris disks in and 
we summarize our results in 

2. DATA ANALYSIS 

Our data comprise monthly timing observations with 
the Parkes telesc ope of 151 pulsars , large ly identical 

The bulk 
the 20 cm 


to those given by Weltevrede et al. 

carried 


( 20101 . 


with 


of these observations are carried ou 
multi-beam receiver, with 256 MHz of bandwidth cen¬ 
tred at 1369 MHz digitally polyphase filterbanked into 
1024 channels and folded in real-time into 1024 phase 
bins. Observation lengths fobs range from 2-20 minutes 
depending on the flux density and sharpness of the pulse 
profile. Data typically span the seven years from 2007 
June to 2014 Oct, and within this interval we supplement 
our observations with archival data where available. 

For each pulsar, we have determined the best timing 
model, typically comprising spin-down terms v and i>, 
the pulsar position, and, where required, proper motion 
and glitch parameters. A representative sample of the 
timing residuals to such a model (including a v term) 
appears in Figure The pulsars are strongly affected 
by red “spin” noise, and to estimate these parameters 
robustly requires modelling the resulting covariance be¬ 
tween TOA measurements. We assume the red noise pro¬ 
cess is widesense stationary and model its power spectral 


P{f) = Ao [1 + U/Uf 


I -a/2 


( 1 ) 


Here, /c allows a possible low-frequency cutoff to the tim¬ 
ing noise, while a describes its shape. By the Wiener- 
Khinchin theorem, P{f) is the Fourier transform of C(t), 
the covariance between TOAs separated in time by an in¬ 
terval r. To this C'(t) we add the white noise model, a 
diagonal covariance matrix comprising measurem ent un¬ 
certainty and “jitter” noise (Shannon et al.|2014), which 


we model as cr| = /3/tobs.j- The complete covariance 
matrix C can be written as a Cholesky decomposition, 
C = LL*, with L a lower triangular matrix. If the noise 
model is accurate, L transforms the TOAs into uni t vari¬ 
ance normal random variables (Coles et al. 20111. The 
timing noise and jitter parameters, together with the pul¬ 
sar timing parameters, form a complete model, and with 
the Cholesky decomposition we can efficiently evaluate 
its log likelihood. We apply Monte Carlo Markov Chain 
methods (MCMC, |Forenian-Mackey et al. 20131 to de¬ 
termine its maximum and shape, yielding estimators for 
the parameters and their statistical uncertainty. We gen¬ 
erally Hnd this approach results in an acceptable fit, as 
evidenced by the successful whitening of the residuals, 
justifying our adoption of a power-law red noise model. 

To characterize the reflex motion due to a single planet, 
we add to the above model the hve Keplerian parameters 
(orbital period Pb, projected semi-major axis a, epoch of 
periastron Tq, eccentricity e, and longitude of ascending 
node H) that fully define the orbit. We restrict the pa¬ 
rameter space to 0 < e < 0.3 (with uniform prior) and 
56 < Pb (d) < 5000. This is roughly the range of eccen¬ 
tricity expected from oligarchical planet formation scat¬ 
tering (see The lower Pb bound is critically sampled 
with our monthly observing cadence, while the upper 
bound is roughly twice the data span. This is a substan¬ 
tial parameter space, and to ensure we explore it fully, 
we divide it into three subspaces: 56 < Pb (d) < 341, 
341 < Pb (d) < 393, and 393 < Pb (d) < 5000, a division 
motivated below. Following a lengthy burn-in period for 
the Markov chains to “forget” their initial conditions, 
we draw 1.2 x 10 ® samples from the first subspace and 
6 X 10® samples from the latter two to form our final 
MCMC sample. 

To search for planets, we computed the change in the 
log likelihood Jlog£ between the null (no planet) model 
and the best-fitting planet model over Pb bins. The value 
Jlog£ for a given planet mass depends strongly on the 
level of white and red noise, but with simulations we es¬ 
tablished 5 log £ « 12 as a reasonable universal threshold 
for a signihcant detection. 

Restricting attention to candidates with Pb <2 yr, we 
found significant modulation in five pulsars,with the tim¬ 
ing residuals of one, PSR J1705—3950, appearing in Fig¬ 
ure There are two strong arguments suggesting these 
are not pulsar planets. First, the implied masses are 
as great as 6 M 0 , and such objects would be easily de¬ 
tected, but have not been, in sa mples including older pul¬ 
sars w ith little timing noise (e.g. Hobbs, Tyne & Kramer] 
2010 1 . Second, three show evidence tor a signihcant sec¬ 
ond harmonic, similar to the well-known state switching 


PSR B1828—11 (Stairs, Lyne fc Shemar|2000 Lyne et al. 


20101 , inviting a magnetospheric interpretation. It these 
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Figure 3. The 90% containment interval of the posterior dis¬ 
tribution for companion mass, for five representative pulsars, as a 
function of companion period. Symbols indicate the masses and or¬ 
bital period for inner solar system planets Mercury, Venus, Earth, 
and Mars. 

harmonics are instead interpreted as two-planet systems, 
the resulting 2:1 mean motion resonance is dynamically 
unstable for the m ass ratios (Beauge, Ferraz-Mello & 
Michtchenko 2003). We will present a detailed analy¬ 
sis and interpretation of these modulations in a future 
work, but emphasize here the key conclusion, that these 
modulations are unlikely to represent bona fide planets, 
and that the small number of candidates ensures that the 
upper limits we compute below remain valid. 

3. UPPER LIMITS ON PLANET MASSES 

The Keplerian parameter a is related to the companion 
mass by the mass function: 


iTic smz 
Mm 


1.50 X 10-3Its 


M, lyr\ 

Mq pJ 


2/3 


( 2 ) 


where we have assumed iric M*, with M* the neu¬ 
tron star mass. While masses of isolated neutron stars 
are generally unknown, measurements from unrecycled 
binary systems suggest the Ch andrasekhar ma ss M* = 
1.4 Mq is representative (e.g. Lattimer 2012). If the 
plane of the putative binary is randomly orientated, i.e. 
p(i) = sin/, then the mean and median values of sin / are 
7r/4 = 0.785 and 0.866 respectively. We report projected 
mass values rric sin/, so physical masses corresponding 
to our upper limits will typically be 15% greater. 

Several examples of the resulting planet limits, as a 


function of Pf,, are shown in Figure 


^ As noted by 
Thorsett & Phillips| (|1992|) , the generalfeatures of these 
spectra reflect simple physics and the relative strength 
of the white and red noise. The high-frequency limit is 
set by white noise in the data, in turn set by the observ¬ 
ing cadence and the TOA precision. In the absence of 
red noise, this limit scales as following Equation 

The sinusoid induced by planets with Pf, « lyr can 
oe absorbed by shifting the apparent pulsar position, re¬ 
sulting in a dramatically reduced sensitivity. However, 
because we generalize our model to include 0 < e < 0.3, 
sensitivity at these time-scales is actually improved rel- 


Figure 4. Limits on planet masses for pulsars with > 0.3 Myr 
(solid lines) and 0.03 < Tc / Myr < 0.3 (dashed, faint lines). The 
lines indicate the fraction of the posterior distribution of compan¬ 
ion mass, for each period bin, contained below. The grey points 
correspond to confirmed exoplanets, pri marily from Keple r, and 
are taken from the Exoplanet Database i jWright et al.|2011[ |. 

ative to the circular orbit case. Finally, the partial si¬ 
nusoid of a long-period planet is similar to a parabola 
and is degenerate with v, and the realization of tim¬ 
ing noise, admitting large values of a. Our division into 
three subspaces reflects this behaviour. 

Figure shows the range of typical cases. 
PSRs J1028-5819 and J1614-5048 are both af¬ 
fected by red noise and their sensitivity decreases 
monotonically with increasing Pb, but the low white 
noise level for J1028—5819 allows good sensitivity to 
short-period objects like Mercury. PSR J1016—5819 
has a relatively high white noise level but very little 
red noise, and the favorable scaling of sensitivity with 
Pb yields good sensitivity to objects with Pb > lyr, 
like Mars. Finally, J1845—0743 offers low noise at all 
frequencies and provides good sensitivity to every inner 
solar system analogue. 

We can extract a posterior probability distribution for 
the entire population by assuming each pulsar is equally 
likely to host a planet. Then the posterior distribution 
of companion mass in a period bin is simply the sum 
of the individual distributions. Because we do in fact 
expect different pulsar ages to correspond to different 
stages of disk activity and planet formation, we have di¬ 
vided our sample into three bins of characteristic age {tc 
in Myr): Tc > 0.3 (47 pulsars); 0.03 < Tc < 0.3 (84); and 
Tc < 0.03 (20). Our primary result appears in Figure 

t where we have constructed posterior distributions for 
e first two Tc bins. The curves show the companion 
mass below which the indicated fraction of the poste¬ 
rior distribution is contained. They can be interpreted 
as limits on the occurrence of planets of a given mass: 
there is only a 5% chance a planet as massive as those 
indicated by the uppermost curve could reside in our 
sample, etc. For completeness, we performed a similar 
analysis using a frequentist method (profile likelihood) 
but assuming circular orbits, and we obtained similar re¬ 
sults when considering the 95% confidence upper limits 
of the population. 
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The oldest pulsars provide the tighest constraints, and 
these are the most likely to harbor planets, having had 
sufficient time for disk evolution and planet formation. In 
this population, we conclude that planets of O.I-O.2M0 
and Pb <lyr are absent from at least two thirds of the 
sample, and planets with M> 0.6 Mq can occur in no 
more than 5% of systems. Consequently, planet forma¬ 
tion within «1.4AU of neutron stars is at best a rare 
phenomenon. 

Younger systems (0.03 < < 0.3) provide similar but 

somewhat poorer constraints: planets exceeding 0.3 M0 
are absent from at least 2/3 of our sample, and planets of 
mass 1-1.5 M0 must be rare, present in no more than 5% 
of cases. The youngest systems (tc < 0.03) have low sen¬ 
sitivity due to their substantial timing noise. We expect 
such systems are also too young for planet formation (see 

S , but for completeness, we note that our limits rule out 
nets of 1-2 M0 in 68% of the very young sample. 

4. IMPLICATIONS FOR DEBRIS DISKS 


Our upper limits clearly show that Earth-mass planets 
are rare or absent within the first ^2 Myr of the forma¬ 
tion of a pulsar. What implication, if any, does this re¬ 
sult have for the presence of debris d i sks ar oun d neutron 


stars? Fol l owing Phinney & Hansen (1993) and Miller & 


Hamilton (2001), we consider it unlikely that any original 


planets with periastron < 1AU could have survived im¬ 
mersion in the stellar envelope or remained bound after 
the supernova kick, and we concentrate on in situ planet 
formation. 

First, we must ask: does our sample trace the young 
pulsar distribution? Our sample excludes magnetars and 
other obviously young but low-A neutron st ars, e.g. the 


weakly magne tized comp act central objects (Halpe rn & 
Gotthelf|2010 ). However, Watters & Romani (2011) and 


other authors have argued that 7-ray pulsars, with prop- 
erties similar to our sample, require the bulk of Milky 
Way supernova activity to account for their numbers. We 
therefore believe our selection of high-A pulsars presents 
a minimally biased view of the young pulsar population, 
but it does not constrain disk formation around less typ¬ 
ical neutron stars. 

Secondly: given a disk with sufficient mass and an¬ 
gular momentum (radius), does planet formation occur 
efficiently and produce palpable planets? Fallback disks 
differ from protoplanetary disks in a key aspect: there 
is no external source of material. Thus, we appeal to 
a study aimed at reprod ucing the properties of planets 
around PSR B1257-I-12. |Hansen, Shih & Currie (20091 


found that disks with fallback-like angular momenta and 
layered viscous accretion can produce Earth-mass plan¬ 
ets rapidly (10^-10^ yr) at inner solar system scales and 
with relatively small (e < 0.2) eccentricities. Even more 
encouragingly, they found that annular density profiles 
produced planets of a few M0 in < 1 Myr. Such profiles 
might result if dust agglomerates into grains quickly af¬ 
ter condensation during disk expansion, in which case 
the disk can deposit most of its dust at the same ra¬ 
dius. In these scenarios, two or three planets with mass 
1-3 M0 populate the inner AU, along with a few resid¬ 
ual planitesimals, and would be easily detectable in our 
sample. 

Thus, we can say that if such disks are common around 
young neutron stars, at least some should host Earth- 


mass planets and we should detect such in our sam¬ 
ple. Moreover, even if rapid sedimentation does not pro- 
ceed and the mass profile re mains more homogeneous. 


Hansen, Shih & Currie (20091 find the production of de- 
tectable planets is still inevitable and merely proceeds on 
a 1-10 Myr time-scale. The absence of Earth-mass plan¬ 
ets in both our sample a nd in larger samples includin: 


2010 


pulsars with > 10^yr (Hobbs, Lyne & Kramer 
seems to rule out such disxsT 
However, such disks seem to play an important role, 
at least in the pulsar literature. Debris disks have 


ration (Cordes & Shannon|2008 

Brook et al. 12014 

Kerr 

et al. 2 u 14|), altering the spin-d 

own 

histories ot pulsars 

(lAlenou, Perna & Hernquist 

2001a 

Yan, Perna & So- 

ria 21)121, contributing to pulsar timing noise (|Shan- 


non et ^2013|), and even enabling the pulsar iiiecha- 


nism (|Mi chel y Desslei'll 1981 1 ). 


Narayan 


among others. 


Chatterjee, Hernquist & 
lave speculated that nor¬ 


mal pulsars accreting from a fallback disk may provide 
the X-ray emission observed in anomalous X-ray pulsars. 
They may also be responsible fo r some cases of peculia r 
supernova light curve evolution (|Dexter & Kasen||2013 ). 

Can we save this picture by forming disks that evade 
our upper limits? After an initial period of super- 
Eddington accretion, a fallback disk evolves self-similarly 
provided it remains viscous, transporting angular mo¬ 
mentum and a small fraction of the material beyond the 
tidal disruption radius where planet formation can oc- 


ley 

1991) 

SO 1 

ong as 


19911 is generally thought to provide this viscosity 


ong as the disk remains ionized, but as it cools, ther¬ 
mal ionization instabilities resu lt in rapid condensation, 


leaving an i nactive dust disk (Menou, Perna & Hern¬ 
quist 2001b| ). (Though we note accretion may continue 


via an a ctive central “layer”, in which case Currie 


Hansen (2007) find the disk can reach the right scale 
for planet tormation.) The final radius of such a de¬ 
bris disk depends sensitively on metallicity, and high- 
metallicity disks generally became inactive at small radii 
(<10^°cm), within the tidal disruption radius. Such 
disks could still form planets, but would first require 
scattering to move material beyond the tidal disruption 
radius. Moreover, this process is impeded by ablation 


of small bodies by the energetic pulsar wind (Miller & 


Hamilton 2001) and cannot proceed until the pulsar’s 


spin-down power E drops low enough. 

On the other hand, near impact on a companion from 
a supernova kick c ould form a debris disk w ith large an¬ 
gular momentum (Phinney & Hansen 1993), giving rise 
to planets with periods too long to be rOiably detected in 
our analysis. We expect such events to be rare, as even in 
tight binaries the solid angle presented by an unevolved 
companion is small, and simulations suggest most bina - 
ries will have Pb > lOd ( Terman, Taam fc Savage|1998[ ). 
Our sample contains only two binary systems, support- 
ing the idea that most sample members did not form in 
close binaries. 

Finally, it may be that long-lived debris disks simply 
fail to form. Stellar material may fall back onto the 
compact rem nant but ultimately be ejected b y a pro¬ 
peller phase (Ek§i, Hernquist & Narayan! 2005). Or, the 

tha 


pre-supernova conditions may be such that the specific 
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angular momentum of the stellar core is too low to circu¬ 
larize. In this case, the explosion energy either unbinds 
the entirety of the star, leaving a neutron star, or a huge 
amount of material falls back, triggering a collapse to 
a black hole. Recent one-dimensional numerical simula¬ 


tions by Perna et al. (2014) suggest that magnetic cou¬ 
pling saps the specihc angular momentum of the stellar 
core to such a degree that only fine-tuned explosions can 
produce appreciable fallback onto a neutron star, and the 
resulting disk lasts only a few hundred seconds. 

5. SUMMARY AND CONCLUSION 

We have searched a large sample of young pulsars for 
periodic modulation characteristic to planetary compan- 
ions. Our work is an imp rovement on previous efforts 


(Thorsett & Phillips 1992), as our pulsar sample is two 
orders of magnitude larger and we employ sophisticated 
methods to mitigate pulsar timing noise and model re¬ 
alistic (noncircular) orbits. Despite the good sensitiv¬ 
ity to low-mass planets we find no compelling evidence 
for such systems. We argue that such companions could 
have formed in debris disks within the 2 Myr age range 
spanned by our sample, and their absence implies super¬ 
nova fallback disks are either rare or confined to small 

radii. _ 

Whence, then, t he planets of PSR B1257-I-12? |Miller| 
& Hamilton (2001) have argued PSR B1257-I-12 was born 
with its present 6 ms spin period and an intrinsically 
weak magnetic field. Such a weak field may point to 
atypically weak magnetic coupling in the progenitor s tel- 
lar core, evading the constraints of Perna et al. (2014). If 
this is so, isolated millisecond pulsars with modest spin 
periods may make excellent targets for direct debris disk 
searches. 


The Parkes radio telescope is part of the Australia 
Telescope, which is funded by the Commonwealth Gov¬ 
ernment for operation as a National Facility managed by 
CSIRO. 
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